Investigation and comparison of multi-state and two-state atom laser output-couplers 
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We investigate the spatial structure and temporal dynamics created in a Bose-Einstein conden- 
sate (BEC) by radio-frequency (RF) atom laser output-couplers using a one-dimensional mean-field 
model. We compare the behavior of a 'pure' two-state atom laser to the multi-level systems demon- 
strated in laboratories. In particular, we investigate the peak homogeneous output flux, classical 
fluctuations in the beam and the onset of a bound state which shuts down the atom laser output. 
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INTRODUCTION 

The analogy between atom lasers and optical lasers is 

n 

strong 1]. Both optical and atom lasers create a coherent 
output beam of bosons that are photons in the case of an 
optical laser and de Broglie matter waves in the case of an 
atom laser. The lasing mode in an optical laser is pumped 
through a non-thermal equilibrium process and is not the 
lowest energy mode of the cavity but, rather, is a highly 
excited mode many wavelengths long. For atoms, the 
lasing mode is populated through Bose-Einstein conden- 
sation and is the ground state of the trap. In an optical 
laser, a beam is outcoupled using a grating or a par- 
tially reflecting mirror. Most atom laser beams studied 
to date have been outcoupled either continuously (long 
pulse) or pulsed using RF radiation to transfer atoms 
from a trapped magnetic sub-state to a state that does 
not interact with the trapping field 0, EJ 0] ■ The atoms 
fall away from the trap under gravity producing a coher- 
ent matter wave [HEHH- 

Optical lasers have found broad applications in preci- 
sion measurements to address questions both fundamen- 
tal and applied in nature. In many cases, we expect to be 
able to perform such experiments more effectively and to 
higher precision with atom interferometry |3, |lfj ■ Mea- 
surements made with the low density, highly coherent 
beam from an atom laser may not be limited in precision 
by the mean-field interaction that plagues interferomet- 
ric measurements made with more dense atomic sources 
such as full Bose-Einstein condensates 0, ^3 • 

We have recently found that the RF output-coupler has 
a number of properties which hinder the production of a 
high flux, shot-noise limited atom laser [lilll^l . We have 
found that a trade off must be made between output flux 
and classical (or dynamical) fluctuations in the output 
beam. We have measured that the peak homogeneous 
flux into the atom laser beam is significantly below the 
flux that can be provided by the finite reservoir of BEC 
atoms in a typical experiment. This homogeneous flux 
limit is imposed by the interaction between the inher- 
ent multiple internal Zeeman states of the magnetically 



confined atoms. Furthermore, we have also found that 
a previously predicted effect known as the bound state 
of an atom laser effectively shuts off the RF output- 
coupler and hence the atom laser beam. 

Recently, two-state atom laser systems have been pro- 
duced experimentally 0, llfil | with the aim of avoiding 
classical fluctuations existing in more complicated multi- 
state systems. Y. Le Coq et al. produce a 87 Rb BEC 
in the F = l,rriF = — 1 state in a harmonic trap with a 
large bias magnetic field (40G). An atom laser is created 
by applying an RF coupling resonant with a transition 
from the wif = — 1 state to the mp = state. Transi- 
tions to the anti-trapped mj? = 1 state are suppressed 
by the second order Zeeman shift, creating an effective 
two-state atom laser. A. Ottl et al. also create a 87 Rb 
BEC in the F — l,mp = —1 state. The authors then 
use a microwave transition at 6.8GHz to resonantly cou- 
ple atoms to the F = 2, = state. The magnetic 
trap bias field splits the Zeeman states in the F=2 man- 
ifold by approximately 1MHz, once again leading to a 
two-state system. 

In this paper we utilize a one-dimensional numeri- 
cal model to investigate the spatial dynamics of an RF 
output-coupler for five, three and two-level atom laser 
systems. The five and three-state systems correspond to 
the experimentally relevant Zeeman levels of the F=2 and 
F=l ground states of 87 Rb and the two-state simulations 
to the recent experiments of Y. Le Coq et al. and A. Ottl 
et al. The question we answer is whether the significant 
amount of experimental effort to make a two-state sys- 
tem has an actual influence on the quality of the atom 
laser. We analyze each of these systems with respect to 
experimentally important properties. In particular, we 
investigate the peak homogeneous output flux, classical 
fluctuations in the beam and the onset of a bound state 
which shuts down the atom laser output. 

We find that, in the weak coupling regime, the three- 
state and pure two-state system are indistinguishable 
with respect to classical fluctuations and flux. We con- 
clude that the added experimental complications of pro- 
ducing a truly pure two-state system do not represent a 
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significant benefit. 

THE MODEL 

An RF output-coupler uses the resonance between non- 
degenerate Zeeman states in the small bias field at the 
minimum of the magnetic trap. Applying a monochro- 
matic RF magnetic field resonant with this splitting 
coherently couples atoms between the relevant Zeeman 
states. The magnetically trapped F = 2, mp — 2 atoms 
can then be 'shunted' through the F = 2, mj? = 1 state 
to the F = 2,mp = state in which they no longer 
interact with the magnetic potential. The presence of 
the trapping magnetic field also introduces a spatial res- 
onance associated with a given frequency allowing the 
centre of the resonance to be tuned within and around 
the condensate. The resonance condition is satisfied on 
the surface of an ellipsoid centered around the minimum 
of the magnetic field. Gravity introduces an asymmetry 
such that the minimum of the trapping potential given 
by \mtu 2 p 2 is shifted down vertically by G shift — g/u 2 
from the centre of the magnetic field minimum. Here m 
is the atomic mass, g the acceleration due to gravity and 
ujp the radial trapping frequency of the F = 2,mp = 2 
state. This asymmetry produces a 'preferred' direction 
of the coupling process, and the atoms fall out of the 
trapping region. 

Time dependent equations 

Ballagh et al. |l7J introduced the Gross-Pitaevksii 
(GP) equation as an effective tool for investigating this 
type of atom laser within the semi-classical mean-field ap- 
proximation. A number of groups have fou nd g ood agree- 
ment between theory and experiment jlSl lid . |20j , using 
mean-field models of the atom laser. Numerical descrip- 
tions of an atom laser within the mean-field frame work 
are complicated by the large velocities that atoms reach 
when falling in a gravitational potential. The resultant 
small de Broglie wavelengths require very fine temporal 
and spatial numerical grids in order to accurately fol- 
low the dynamics. Additionally, to run simulations that 
reflect experiments on a time scale longer than a few mil- 
liseconds, an apodising mechanism must be introduced to 
absorb the beam and hence avoid breakdown of the nu- 
merical techniques used. These factors make simulating 
and understanding atom laser dynamics a complicated 
proposition. In order to simplify the numerics, a true 3D 
system can be transformed to lower dimensions. The di- 
mensionality reduction can be performed non rigorously 
by writing an equivalent equation for the system in the 
dimension(s) of interest The F=2 GP model 0] 

of the atom laser in one dimension (the coordinate over 
which gravity acts) is then given by 



i<p 2 = (£ + z 2 + Gz- 2A)0 2 + 20^i 

ifa = {C + ^z 2 + Gz- A)0i + 20^2 + V6ft0 o 

i4>u = {c + Gz)ct> + Vficifa + v^fty-i 

i4>-! = {C- ^z 2 + Gz + A)0_i + 2f20_ 2 + V6fl(j>o 
= {C-z 2 + Gz + 2A)0_ 2 + 2f^! , 

(1) 

where <pi is the GP function for the zth Zeeman 

state and C = -§ + U (Sf = _ 2 \^i\ 2 ) . Here A and 

are respectively the detuning of the RF field from 
the Bq resonance, and the Rabi frequency, measured 
in units of the radial trapping frequency uj z (for the 
F = 2, mp = 1 state), U is the interaction coefficient 
and G = w^-i— — W 2 gravity. The wave functions, time, 
spatial coordinates, and interaction strengths are mea- 
sured in the units of (Ti/mii^) -1 / 4 , w" 1 , (H/muJz) 1 ^ 2 , 
and (hu! z )~ 1 (H/mLij z )^ 1 ^ 2 , respectively. The nonlinear 
interaction strength is derived by requiring that the ID 
Thomas-Fermi chemical potential be equivalent to the 
3D case. There are no free parameters in this model; we 
use U = 6.6 x 1(T 4 , G = 9.24, f2 = - 14, A = 10.7. 
Although all the numerics are executed in dimensionless 
units, data are presented in dimensional units. 

Fortuitously, the gp factors for the F=2 and F=l 
ground states of 87 Rb are the same except for a change of 
sign. The three-state and two-state atom laser systems 
are then simply subsets of the five-state equations pre- 
sented above, with the modification that the pre-factor 
\f& in the coupling terms between the mp — —1,0,1 
states becomes \ 2. 

Our simulations are performed in a commercially avail- 
able and widely used platform Matlab, using a fourier 
based, symmetric split step algorithm. For this work we 
use a 2048 point spatial grid from -40 to 40 with the equi- 
librium position of the condensate at 20 and a temporal 
resolution of 10 -4 . Apodising boundaries for each state 
have an exponential form and are positioned in order that 
no spatial aliasing occurs. 

In addition to solving the time-dependent GP equa- 
tions for the atom laser we also find the time- independent 
solutions to provide accurate initial conditions for our 
simulations. The time-independent GP equation can 
be found from Eq^ using the substitution </> 2 (r, i) = 
tp(Y)e l ^ t . We use a relaxation technique 0] which, if 
required, allows us to find excited stationary states as 
well as the ground state solution. For the five-state sys- 
tem we use as our initial condition the solution to the 
following time-independent equation, 

^ = {-\^ + u 2 \4>2\ 2 + z 2 )4> 2 . (2) 
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For the two and three-state systems we solve 

Mi = (-^ + tfil<M 2 + ^>-i. (3) 

In these equations the one-dimensional interaction 
strengths are calculated by requiring that the one and 
three-dimensional chemical potentials are equal. 



RESULTS 

In this section we present the general results for each 
of the five, three and two-state systems. Unless other- 
wise stated, for each system studied the RF coupling res- 
onance is arranged to be at the center of the trapped 
BEC by selecting an appropriate value for the detuning. 



Five-state system 

We first analyze the five-state system. In figure 1 we 
show the population dynamics for strong and weak cou- 
pling. Although we are using an apodising boundary in 
the numerical simulations we keep track of the total pop- 
ulation in each state by calculating at each timestep the 
one-dimensional flux passing through a point (we like to 
think of it as a detector) on the numerical grid, signifi- 
cantly below the trapped condensate and above the apo- 
dising boundary. Hence, the numerical grid is effectively 
divided into two sections: above and below this 'detec- 
tor'. At any time during the simulation we can determine 
the total number of atoms in a particular state by sum- 
ming the number of atoms on the grid above the detector 
with the total number of atoms that have passed through 
the detector. Apart from giving results that are more in- 
tuitively obvious, this method allows us to monitor the 
normalisation of the numerics. A number of theoretical 
works have suggested that the atom laser may have a 
'bound' eigenstate 0, H3, H3 , based purely on the exis- 
tence of coupling between a single trap mode and a con- 
tinuum of un-trapped states. Furthermore, in the context 
of producing a two-dimensional BEC, it has been shown 
recently that trapping of all mp states is a natural con- 
sequence of combining RF coupling with a DC magnetic 
trap [3 EH- This trapping can be understood by consid- 
ering the 'dressed-state' basis in which the RF coupling 
and DC potentials seen by the atoms are diagonalised. 
In this basis the dressed eigenstates are linear combina- 
tions of the bare Zeeman states, trapped in effective po- 
tentials created by the avoided crossings. Assuming the 
strong coupling limit and a sudden non-adiabatic pro- 
jection onto the dressed-states, diagonalization yields a 
prediction of up to 62.5% of the initial condensate atoms 
remaining trapped for the F=2 atom laser (four of the 
five dressed-states allow some trapping). 



In figure 1(a) the condensate (trapped atoms in the 
mp = 2 and mp = 1 states) decays slowly and mono- 
tonically, with some small modulation introduced by the 
output-coupling process itself. This behavior is indicative 
of the intermediate output-coupling regime in which most 
experiments have been operated In contrast, in 

figure 1(b) we see the formation of a bound state in the 
case of strong output-coupling. At t = the high power 
RF coupling is switched on and the mp = 2 Zeeman state 
is projected onto the new dressed-state basis as discussed 
above. After a short (~ lms) high frequency exchange a 
fraction of un-trapped dressed-states are ejected, leaving 
the remnant dressed-states to oscillate in the magnetic 
trap. 

Although the population dynamics in figure 1(a) is 
smooth, wc anticipated that the details of the spatial 
dynamics would not be. This is because in the five- 
state system atoms trapped in the mp — 2 state must 
pass through the mp — 1 state to get to the un-trapped 
mp = state. Because atoms in the mp = 1 state have 
a different equilibrium position to the mp = 2 state (the 
gravitational sag is different for each state) the mp = 1 
atoms start to oscillate in the magnetic trap. The dynam- 
ics of such an oscillation is shown in figure 2. This phe- 
nomena is independent of the output-coupling strength, 
and indicates that even at low flux a five-state atom laser 
system will be modulated by classical noise. Density fluc- 
tuations in the output beam are shown figure 3, for dif- 
ferent Rabi frequencies in the weak coupling regime. We 
note here that the spatial dynamics imposed on the atom 
laser beam by the mechanism of the mp = 1 sloshing 
in the trap are an independent noise mechanism com- 
pared to the back-coupling dynamics investigated previ- 
ously 0. 

Three and two-state systems 

The three-level system offers the possibility of a cleaner 
output than the five-level, as there is no intermediate 
state between the trapped condensate {mp = — 1 Zee- 
man state) and the un-trapped beam {mp = Zeeman 
state). Fluctuations in the output will then be due to 
the back-coupling and depletion of atoms to the anti- 
trapped mp = 1 Zeeman state as observed previously 
14]. However, in the limit that the output-coupling be- 
comes weak these effects will be negligible and the sys- 
tem should produce a classically quiet atom laser beam. 
In figures 4 and 5 we show the population dynamics for 
the weak and strong coupling regimes for the three and 
two-state systems respectively. The behavior of the three 
and two-state population dynamics is qualitatively simi- 
lar to the five-state system. For weak coupling strengths 
the condensate decays monotonically, with a greater frac- 
tion of condensed atoms transferred to the mp = state 
than the five-state system under similar conditions. In 
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the strong coupling limit the behavior of the system is 
similar to the five-state, although here there is a more 
clean cyclic oscillation of atoms between the three Zee- 
man states. Once again a percentage of all states remain 
trapped. 

The oscillatory population dynamics observed in fig- 
ures 1, 4 and 5 (b) for the strong coupling limit is also 
reflected, if not driven by, periodic spatial oscillations. 
In figure 6 we show an example of such oscillations for 
the five-state system. In this figure one see the ejection 
of the un-trapped dressed-states early in the simulation, 
and then a clean periodic oscillation. Two oscillation pe- 
riods are clear from this figure, corresponding to the F=2 
and F=l radial trapping frequencies. Interestingly, these 
oscillations have their upper maximum positions at the 
minimum of the magnetic trapping potential (the point of 
gravitational sag for the F=2 state), and hence oscillate 
up only one side of the total potential (trap plus grav- 
ity) . Similar behaviour occurs for the two and three-state 
systems as well. We believe that further detailed theo- 
retical and experimental study of these unusual dressed- 
state spinor (multicomponcnt) condensates will reveal a 
rich and potentially useful phase space. Already, dressed 
systems have been used to demonstrate two-dimensional 
trapping potentials |2^| and to investigate double well 
interferometry with condensates |2( : il| . 

Comparison between two and three-state systems 

In this section we compare the dynamics of the two 
and three-state atom lasers. As mentioned in the intro- 
duction, significant effort is required in order to produce 
a 'closed' two-state atom laser as the relevant alkali atom 
manifolds have at least 3 Zeeman states linked by allowed 
RF transitions. Given that back-coupling fluctuations 
and the bound state arise even in the two-level system, 
limiting the homogeneous output flux as in the 'natural' 
systems, it is prudent to ask whether the two and three- 
state atom lasers actually differ that much in the weak 
coupling limit. 

In figure 7(a) we compare the time dependent density 
of the two and three-state atom lasers at a point below 
the condensate. It is clear that, for the one-dimensional 
model we are using, the two systems are essentially indis- 
tinguishable in the weak coupling regime. It is interesting 
to note that even for stronger coupling, the details of the 
classical fluctuations on the beam are also mirrored in 
the two systems. The difference in amplitude between 
the simulations is accounted for by the loss of atoms to 
the anti-trapped Zeeman state in the three-level system. 
In the weak coupling limit these atoms are expelled from 
the system on a time scale much faster than the back- 
coupling time set by the Rabi frequency, and so have 
little effect on the dynamics of the system. 

In figure 7(b) we show the density fluctuations as a 



function of time for the two or three-state atom laser. 
It is clear that even in the weak coupling limit taken 
here the density fluctuations due to back-coupling are 
a significant contribution to the noise in the atom laser 
beam. This has major implications for the use of atom 
lasers in precision measurement systems, where high flux 
and minimal classical noise are essential features required 
of the atom beam. 

In figure 8 we compare the number of atoms in the 
condensate to the number of atoms in the mp = atom 
laser state for 15 ms of output-coupling. There is a clear 
peak in the number of atoms transferred into the mp = 
state as a function of Rabi frequency for all systems con- 
sidered in this paper, with the two-state system having 
the highest coupling efficiency and the five-state having 
the lowest. However, a comparison of figure 8 with fig- 
ure 7 reveals that we will not be able to operate the 
atom laser near the peak output-coupling rate because 
of increasingly severe density fluctuations, which impede 
the usefulness of the beam for measurement applications. 
For example, the peak in the two-state system is around 
il = 500Hz while the density fluctuations are already 
significant at a Rabi frequency of 130Hz. 

In summary, we find that the peak homogeneous 
output-coupling rate achievable in an atom laser to be 
significantly lower than the maximum output-coupling 
rate. For the one-dimensional model considered here we 
find that, in this homogeneous output-coupling regime, 
there is practically no difference in our results between 
the three and two-state atom lasers. One may ask 
whether the small amount of anti-trapped atoms gen- 
erated in the three-state system will effect the results of 
a correlation function measurement on the atomic beam, 
such as the one made recently in the Zurich group |l6| . 
In this experiment an atom laser beam is produced in 
the weak coupling regime and subsequently falls through 
a high finesse cavity, in which the researchers are able 
to detect the passage of a single atom with high tempo- 
ral resolution. A simple calculation shows that the small 
amount of anti-trapped atoms produced by a three-state 
atom laser system would miss the cavity by a large mar- 
gin, thus having no effect on the results of the measure- 
ment statistics for the mp = state. One can imagine 
many experiments where a true two-level system is im- 
portant, but we conclude that the atom laser is not one 
of them, at least in the weak coupling limit. 

Discussion of results with respect to future 
experiments 

We plan to lead further experiments in the near future. 
By creating a very large condensate we should be able to 
have both resonance positions (of the mp = 2 and mp = 
1 trapped states) inside the condensate. One can expect 
the spatial dynamics imposed on the atom laser beam 
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by the mechanism of the wif = 1 sloshing would then 
disappear. However, preliminary simulations carried in 
this configuration seem to show no difference and further 
investigation is necessary. 

A Raman output-coupling setup will also be imple- 
mented as it overcomes many problems with RF output- 
couplers [j^. In particular, in the case of a Raman atom 
laser, we could selectively outcouple atoms in a given Zee- 
man state, taking advantage of the spatial dependance of 
the resonances. 

CONCLUSION 

In this paper, we have characterized how RF output- 
coupling affects the spatial structure and temporal dy- 
namics in a BEC. We also presented the ID mean-field 
model used for this purpose. The peak homogeneous out- 
put flux, classical fluctuations in the atom laser beam and 
the onset of a bound state were investigated in the cases 
of multi-state and two-state atom laser systems. The re- 
sults showed the five-state system was clearly inappropri- 
ate to create a clean and homogeneous atom laser beam. 
It also appeared the two and three-state systems were al- 
most indistinguishable in the weak coupling regime where 
the atomic beam is homogeneous, with no density fluctu- 
ations. As a conclusion we think the 'natural' three-state 
system should be preferred as it does not require any ad- 
ditional experimental devices. 

We are currently studying a momentum kicked Raman 
coupling based on the same model and it will be inter- 
esting and important for future measurement devices to 
compare the results to those described in the present pa- 
per. 
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FIG. 1: Comparison of the population dynamics in the five-state system for (a) weak and (b) strong coupling. Parameters are 
U = 6.6 x 10" 4 , G = 9.24, A = 10.7. 
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FIG. 2: Spatial profiles of the five-state system, showing the effect of the mr = 1 state oscillations on the output mode. The 
color axes are adjusted to give the best contrast for each of the Zeeman states. Each plot shows the atomic density, with 
the complete spatial grid in the vertical direction and the temporal grid in the horizontal, (a) Total density (the sum of the 
densities in all five Zeeman states) as a function of time, (b) density of the tuf = 2 state, (c) density of the rriF = 1 state, and 
(d) rriF — state. The dynamics were obtained after 15 ms of intermediate output-coupling (fi = 1kHz). 
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FIG. 3: Density fluctuations in the rriF = state of a five-state system at a single point in the beam as a function of time. 
The fluctuations due to the sloshing of the mp = 1 are increasingly severe as the coupling strength is increased. Parameters 
are U = 6.6 x 10" 4 , G = 9.24, A = 10.7. 
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FIG. 4: Population dynamics in the three-state system for (a) weak and (b) strong coupling. Parameters are U = 6.6 x 10 4 , 
G = 9.24, A = 43. 
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FIG. 5: Comparison of the population dynamics in the two-state system for (a) weak and (b) strong coupling. Parameters are 
U = 6.6 x 10~ 4 , G = 9.24, A = 43. 
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FIG. 6: Periodic spatial oscillations in the population dynamics of a five-state system. One can see the ejection of the un- 
trapped dressed-states as well as both oscillation periods corresponding to the F=2 and F=l radial trapping frequencies. Note 
that the upper edge of the oscillations is at the minimum of the magnetic trapping potential. 




FIG. 7: Comparison of the time-dependent density of the two and three-state atom lasers at a point below the condensate after 
15 ms of output-coupling and for different Rabi frequencies (a), (b) shows the time-dependent density of the two-state atom 
laser at very low intensities after 100 ms of output-coupling. 
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FIG. 8: Comparison of the number of atoms in the condensate and in the rriF = atom laser state for the two-state and 
multi-state systems after 15 ms of output-coupling. 



